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ABSTRACT 



This thesis addresses the problem of determining the 
optimal number of spares for a multi- item inventory system 
with a procurement budget constraint. Various inventory 
models are considered with objective functions like time- 
weighted units short, units short, essentiality -weighted 
units short and pseudo-availability. Solution algorithms 
are derived using the generalized Lagrange multiplier 
approach and a marginal analysis approach. 

Sample data and output results are provided and compari- 
sons of the alternative models are given. Finally, a dis- 
cussion and example is given -of the use of the models as 
a means of estimating the budget required to attain a speci- 
fied level of performance. 
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I. 



INTRODUCTION 



In today's world, while all systems are becoming more 
and more sophisticated, the control and maintenance of 
inventories of these systems is a problem common to all 
enterprises and military services. In private and commer- 
cial concerns the effective control of inventories can re- 
sult in decreased costs, increased sales and profits and 
consumer satisfaction. In the military proper management of 
inventories may contribute to increased availability and 
readiness, decreased inventory investment and system costs. 

For each component of each weapon system two fundamental 
questions must be answered: 

(1) When to replenish the inventory; 

(2) How much to buy for the replenishment. 

In order to answer these questions, many inventory models 
have been developed in the past 30 years. See for example, 
Hadley and Whitin [Ref. 1] , Muckstadt [Ref. 2] and Eriksson 
[Ref. 3] . Most previous work solves a variety of cost 
minimization problems considering expected values of steady 
state variable costs associated with shortage cost, ordering 
cost and storage cost. 

Such models may be appropriate for the commercial sector, 
but are not always appropriate in the military world. 
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In the commercial sector, the objective function of the 
inventory model is to maximize profit or minimize the average 
annual costs. Non-cost oriented objective functions frequently 
are used in the military inventory systems. For example, 
attempts are often made to maximize availability or fill 
rate, or minimize the number of backorders or expected time 
weighted stockouts, or minimize the probability of a stock- 
out with a budget constraint. 

Obviously costs are important in every inventory model. 
However, many real-world inventory problems are so compli- 
cated, one cannot represent accurately the real situation. 

Thus, some simplifications and approximations are used when 
constructing a mathematical model of any real world system. 

If this is not done, the results obtained by use of the model 
can easily lead to operating rules which are worse than those 
currently in use, worse than those which could be derived 
from simple heuristic intuitive considerations. 

Many of the inventory problems are viewed as single period 
problems. For example, initial provisioning, allowance list 
determination and the fly-away kit problem are single period 
problems. These models are perhaps the simplest of the models 
in which demand is treated as a stochastic variable. 

Reasonable objective functions in these models are to 
maximize performance subject to a constraint on the resources. 
Typical measures of performance might be availability, time- 
weighted units short, fill rate, the number of backorders, 
and mean supply response time. 
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This thesis considers various single period models which 
attempt to maximize performance subject to budget constraint. 

Chapter II describes the general single period problem 
and introduces the method used in this thesis of solving 
those problems. 

Chapters III and IV develop the -time-weighted units short 
model and the availability model, and explain the solution 
procedure. Sample data runs for both models are provided. 

Chapter V provides a comparison of models considered in 
the thesis and discusses some of the properties of each 
model, and Chapter VI discusses the use of models for purposes 
of determining the budget required. 

Chapter VII summarizes the results of the research and 
concludes with some suggestions for additional research. 
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II. THE GENERAL PROBLEM 



In this chapter, we consider the general single period 
model as a process for transforming resources into new dis- 
tributions of inventory positions over the line items in 
the inventory . 

The essential problems of control in a line item-inventory 
control system with multiple line items are; 

(1) How much resources to commit at a point in time; 

(2) How shall these resources be allocated to achieve 
system objectives. 

In a typical continuous review inventory system, we can 
determine the optimal order quantity (Q) and reorder point 
(r) for a given item by minimizing the average annual variable 
costs. But in applying this theory to the real world inven- 
tory systems which consist of multiple line items, it is 
frequently the case that resulting minimum cost solutions are 
not feasible because of a budget limitation or some other 
constraint. Thus, in a constrained multi-item inventory 
system, the typical continuous-review policy is sometimes 
inappropriate. In the following section we discuss several 
objective functions to guide the line item inventory control 
system in determining how to allocate availabile procurement 
funds at a particular replenishment epoch. 
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A. GENERAL FORM OF OBJECTIVES WITH CONSTRAINTS 



Consider the case in which an administrator, responsible 
for the replenishment decisions, determines replenishment of 
stocks of various line items on a periodic basis. Suppose 
that a fixed amount of procurement budget has been allocated 
to the replenishment epoch at hand and that a target number 
of reorder actions has been established as a working con- 
straint for the allocation epoch. The administrator's task 
is to transform the available resources into replenishment 
orders for different items. 

1. ' Measure of Effectiveness and Objective Function 

Daeschner [Ref. 5] examined the constrained line-item 
allocation problems. He considered several possible objective 
functions which can be adapted to the case where unsatisfied 
demands are backordered and to the case where unsatisfied 
demands are lost sales. Let tt ^ > 0 be the penalty (reward) 
per unit for item j and let be the demand for item j in a 



period. 


Let 


Xj be 


the 


inventory 


position for 


item j 


after 


ordering 


in a 


period. 


Let D . = i 

3 


d . . Then the 

3 


number 


of sales 


for item 


j in 


the 


period is given by 










d. 


if 


d . < 


X . 
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3 - 












X. 


if 


d. > 


X. 
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The expected 


sales 


for 


item j is 


therefore 
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X . 

f d . p (D . = d . ) + X p (D . > X ) 
d.=l 3 3 3 3 3 3 

3 

which is equivalent to 

oo 

E(D ) - l (d - X )p(D =d.) . 

3 d . =X . +1 33 33 

3 3 

We assume that the inventory system seeks to minimize the 
expected penalty incurred, or, equivalently, to maximize the 
expected penalty avoided. Mathematically, the objective is 
to maximize 



N 



Z(X) = l tt . (E (D . ) - 

j=l 3 3 



l (d -X )p(D =d )) 

d.=X.+l 33 33 

3 3 



Several interpretations and uses of the penalty coefficient 
tt j are possible. Four are illustrated in Table I below. Each 
reflects a formulation of system objectives which has been 
adopted or considered by the Navy Supply System. Daeschner 

[Ref. 5] also considered tt . as a linear combination of various 

3 

coefficients in his line item allocation model. 

There are many other types of objective functions which 
are currently used in the military. For example, 

(1) Minimize units short in a given period 



N 

Z(X) = l l (d - X ) p ( D . =d,) 

j=l d.=X.+l 3 3 33 
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TABLE I 



INTERPRETATIONS 


AND USES OF TT . 

3 


Penalty Coefficients 


Objective 


IT . = c . 

J J 


Maximize expected sales from 
stock. 


tt . = i/y . 
3 3 


Maximize the expected requi- 
sitions filled (y^ = average 
quantity of item j demanded 
per requisition) . 


IT . = 1 

3 


Maximize the expected number 
of units issued from stock. 


7 r . = LT . + TMNIS 
3 3 

- TMISS 


Maximize expected customer 

waiting time per unit avoided 

by issue from stock, where 

LT . is the lead time for item 
3 

j , TMNIS is the calendar time 
anticipated to process a 
request and TMISS is the time 
to affect issue from stock of 
a demand, available item. 



(2) Minimize time-weighted units short 
N 

Z (X) = T TWUS • (X. ) • E. 

j£i 

where: 

TWUS(X^) = time weighted units short 
E^ = essentiality 
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(3) Maximize system availability 



N 

z (x) = n a (x.) 

i=l 1 1 



where : 



= Availability for each item. 

The objective functions (2) and (3) will be explained in 
Chapters III and IV. 

2 . The Line-Item Allocation Model and Solution Procedure 
In the previous sections, many kinds of objective 
functions are introduced. If we define them correctly, those 
objective functions can be solved by various techniques. It 
is evident that an actual inventory system with limited re- 
sources might be unable to carry out a prescribed inventory 
policy if either the amount of procurement funds available 
or the number of replenishment actions exceed the available 
resources. The problem is made more complicated by the fact 
that the objective functions are "non-linear" and the require 
ment that the X^’s must be integers. The problem is stated 
mathematically as 

max Z (s) 

N 

(Al) s . t. £ c . s . < B 

j=i 3 3 

N 

y H ( s . ) < R 

j-1 D ■ 
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s 




f , 



N) : integers 






1 if s . > 0 
3 



= 0 otherwise 



B = the procurement budget limit at the realloca- 
tion epoch 



R = the maximum number of individual procurement 
activities allowed in the present allocation. 



To solve the problem (Al) , the generalized Lagrange multiplier 
(GLM) method of Everett [Ref. 4] can be used. Using this 
method, the problem can be reexpressed as 



£ e s and A^, ^ 0 with optimal solution £* ( A_) . 



Problem (A2) is the Lagrangian problem associated with (Al) . 
Using Everett's theorem, one can determine a bound on the 
optimal solution, Z(s_*) to be 

Z(s*) < Z(s*(A)) - A 1 (B(A) -B) - A 2 (R(A) - R) 



(A2) 



max L(S,AJ 
s 



N 

Z(s) - A. ( ( l c.S.) - B) 
1 j=i 3 J 



N 

- A„( ( l H ( s . ) ) - R) 
j = l 3 
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where 



N 



B ( A) = l c S*(A) 
j=l 3 3 



N 

R(A) = l H(S*(A)) 
j=l 3 



and s*(A) is the optimal solution vector for a given pair 
(A^/A 2 )» In solving problem (A2) , we can separate the N- 
variable optimization problem into N one-variable problems. 
Choosing trial values of A^ and A 2 we maximize 



When considering the integer nature of decision variables, 
the optimal solutions for (A3) are determined by finding the 



(A3) 





values s* such that 
3 





ALj (Sj ,A) 



L. (s . +1 , A) - L. (s . ,A) < 0 

3 3 - 3 3 “ 



In order to get an optimal solution, Daeshner [Ref. 5 ] used 
an interactive computer program, which evaluates the current 



optimal solutions with A^ and A 2 * Each time, the user can 
select a pair A^, X 2 and objective function type to be con- 
sidered. Then the user is provided with output which indi- 
cates the budget consumed, the number of stock replenishments 
generated, the achieved objective function value and a maximum 
attainable value for the objective function. 

After examining the output, the user can modify the 
input parameters and continue or terminate the run. Decreas- 
ing the non-negative multiplier values tends to use more of 
the corresponding resources, increasing the values used, less. 
When the replenishment actions generated by a pair of values 
( X ^ , X 2 ) exactly consume the available resources, B and R, 
the solution is optimal. Frequently exact equality may be 
impossible because of integer nature of the problem. Thus 
the solution obtained may not be optimal, but the difference 
is not likely to be significant. 

B. AUTOMATING SEARCH ON THE LAGRANGE MULTIPLIER 

The interactive search method cannot guarantee an optimal 
solution and it requires trial and error to get the approxi- 
mate optimal solution. Consider the case in which there is 
only a single constraint, with the same type of objective 
functions. The mathematical program is then 

max Z (sO 
s 

(B1) N 

s . t. y c . s . < b 

1-1 1 1 - 
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£ = ( s i , s 2 / • • • ,s N ) = integer number of buys 

B = budget limit 

= price of item i. 

We can rewrite the above equation using a Lagrange multiplier, 
as : 



N 

(1) L(s,0) = Z (s) - 0[ £ c.S. - B] 

j-i 1 1 



Then separate the equation. 



N 

(2) L(s 1/ s 2 , • • • f s N ) = l (Z(s i ) - 0c i s i ) + 0B 



Equation (2) can be maximized by maximizing each sub- 
objective function. If Z(s^) is differentiable with respect 
to each s^, the optimal solution is obtained by 




dZ(s.) 

~ds~. 

i 



0c . , 

l 



i = 1 , 2 , . . . ,N 



Thus set 




= 0 and get 



dz (s ) 

e = ^ — — /c i 



ds 



N 

where 0 is such that Y c.s. = B. Everett [Ref. 4] shows 

j=l 1 1 

that 0 can also be interpreted as a shadow price for the 
objective function: i.e., 0 = 8Z/3B. Due to the integer 
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nature of s^, it is often impossible to get an exact optimal 
solution. Difference equations must be used because the 
region in which the solution is desired consists of a set 
of discrete points. Therefore, let 

(3) AL . (s . , 6 ) = L . (s . , 6 ) - L . ( s . -1 , 0 ) 

l i i i i i 

= Z.(s.) - 0c. s. - Z.(s.-l) + 9c.(s.-l) 

ii ii ii ii 

= AZ^(s^) - c^0 

We know that Equation (3) is a concave function at the point 
£ >_ 0 . The optimal solution must satisfy AL^(s^,0) >_ 0 and 
ALfCSi+l,©) < 0. Thus the optimal solutions are given by 
finding the largest s^'s such that 

AL i (s i ,9) _> 0 



or equivalently 

C4) AZ^(s^) - c^9 £ 0 i = 1 , . . . , N . 

The Lagrangian multiplier 0 can be found by the following 
search algorithm. 
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STEP 1. Find an initial upper bound 0 . Let all s. be 
, u 1 

assigned zero at the beginning and find the 
change of objective function per unit dollar as 
a result of increasing to one unit. 

AZ 1 (1) 

C 1 

AZ 2 (1) 




where AZ^(l) = Z^(l) - Z^(0). 

Because of decreasing marginal returns or objective function 

values and because of the interpretation of 9 , an upper bound 

on 0 is given by : 0 = max [0.,0 9 ,...,0 ] . 

u -L z n 

STEP 2. The initial 0 will be 

o 



0 

o 



+ e u 
2 



where 9^ = 0. 

Find for each i, the largest s^ so that 



AZ i (s i> 
c . 

l 



9 

o 



and evaluate the objective function and the budget 
required. 
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STEP 3. If the budget 
budget, let 



otherwise 

e + e T 

o o L 



Each time update the S vector, the objective 
function values, the upper bound of the objective, and the 
amount of budget consumed. 

STEP 4. Stopping rule. 

Stop when the used budget is equal to the given 
budget or the difference between the current upper 
bound and the objective function value is less 
than some limit (e) . Otherwise go to step 2, 
and continue until the above conditions are 
satisfied. 

A FORTRAN program for this algorithm is given in Appendix B. 

C. MARGINAL ANALYSIS PROCESS 

The theory of marginal analysis has been used in many 
inventory models when resource constraints are active. In 
an economic sense, AZ^(s^)/c^ can be interpreted as the 
marginal increase in the objective function per dollar spent 



used is greater than the given 



( 9 + e ) 

o u 
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achieved by adding one more unit of stock. It is reasonable 
for an inventory controller who has a scarce resource such as 
a procurement budget to buy an item which gives the maximum 
benefit per dollar spent. 

By using a simple computerized algorithm, the line item 
allocation problem can be solved easily. The first step is 
to set all s^ = 0 and compute 



max 

i 



AZ-^s-j+l) AZ 2 (s 2 +1) 



AZ (s+1) 
n n 



n 



j 



If the maximum is taken on for item j , set s . = 1 and 
deduct the unit price for unit j from the budget. The second 
step is then to recompute AZ^ and then find 

AZ. (s.) AZ. (s.) 
max {max { — - — — } , — — 2 — } 
i 3 

The next unit is assigned to the index j where the maximum 
is taken on, etc. This is continued until adding an additional 
unit exceeds the budget constraint. It should be noted, how- 
ever, that the method described does not insure optimality 
[Ref. 1]. Specifically the method may stop too soon. If the 
item i selected from the marginal analysis has a c^ value 
greater than the remaining budget, the procedure terminates 
even though some other item j may have a c^ value less than 
the remaining budget. An obvious improvement in this area 
could be the inclusion of a subroutine that would select 
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from the remaining items the best one from those having c^’s 
smaller than the remaining budget. A FORTRAN program for 
performing this marginal analysis is provided in Appendix C. 



D. SAMPLE DATA RUNS OF UNITS SHORT MODEL 

A weapon system consists of 10 components. The system 
manager wants *to minimize the number of units short by supply- 
ing spare parts to support the weapon system. Suppose that 
the demand rate, lead time, price and essentiality code for 
each item i are known. The objective function can be expressed 
by 



(Cl) 



minimize Z(s) 



10 00 

l l (d -s i )p(D i =d )E i 
i=l d.=s.+l 1 1 1 1 1 



10 

subject to £ s i c i — B 

i=l 



where : 



E^ = essentiality code 

B = budget limit 

-A . T . d. 

e 1 1 ( A ■ T . ) 3 

e (D i =d i> - <a-ri — 

The approximate solution of (Cl) can be obtained by the marginal 
analysis method. Table II shows the computational results 
for this system with known input data. 
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TABLE II 



THE RESULT OF UNITS SHORT MODEL 



No. 


A . 

i 


Lead 

Time 


Price 

($) 


Essen. 


Alloca- 

tion 


Units 

Short 


1 


1.0 


1.0 


10.2 


1.0 


5.0 


0.0007 


2 


0.1 


1.0 


20.0 


1.0 


1.0 


0.0048 


3 


3.0 


1.0 


100.0 


1.0 


2.0 


1.2489 


4 


25.0 


1.0 


2.0 


3.0 


42.0 


0.0013 


5 


1.0 


1.0 


5.0 


1.0 


5.0 


0.0007 


6 


0.5 


1.0 


5.0 


3.0 


4.0 


0.0002 


7 


10.0 


1.0 


1.0 


1.0 


21.0 


0.0012 


8 


5.0 


1.0 


100.0 


1.0 


4.0 


1.4368 


9 


1.0 


1.0 


50.0 


1.0 


3.0 


0.0233 


10 


2.0 


1.0 


100.0 


1.0 


2.0 


0.5413 



Table II shows several properties of the units short 
model. First of all, more than one unit short in a year 
occurs in the high cost items (items 3 and 8) . Second, low 
demands and low price items are allocated enough. Items 2, 

5 and 6 are allocated more than five times their mean demand. 
Also this model tends to stock more of the high demands and 
low price items . 

Finally, the essentiality weights cause greater alloca- 
tions to be provided to those items with high essentiality 
than would be provided with equal weights . 



25 



Other results include: 



Total objective value 3.26 



Shadow price 



0.001899 



Budget limit 



$1170 



Budget left 



$ 0.0 



Az (s i ) E i 



The shadow price is the last maximum value of — . 

i 

It can be interpreted approximately as the amount of 
decrease in the objective function achieved by adding one 
more dollar. 
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III. TIME WEIGHTED UNITS SHORT MODEL 



A. DESCRIPTION OF MODEL 

In the previous chapter we have discussed various objec- 
tive functions and solution methods for the single period 
inventory problem. In the military services, many measures 
of effectiveness have been used to indicate system performance. 
Among these measures are fill rate, availability, mean supply 
response time, the number of stockouts, and time-weighted 
units short. In this chapter we consider a model which 
minimizes time-weighted- units-short (TWUS) . 

Suppose that a weapon system consists of n components 
and the objective is to allocate a given budget for spare 
parts so as to minimize time-weighted-units short for the 
entire system. Assume that 

(1) procurement lead time and repair lead time are known 
constants . 

C 2 ) demands for each installed unit have a known 
distribution . 

(3) the total amount of procurement budget available to 
spend on all components is fixed. 

(4) the objective is to minimize essentiality weighted 
TWUS. Mathematically, the model can be written 

as 
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min l TWUS Is )E 

S 1 = 1 
s . t . 

where : 

TWUS . (s . ) = 

i l 

SLT 

E. 

1 

B 
C. 

l 

In the above problem, if the TWUS is properly defined, this 
model will be solved easily by using the methods explained 
in Chapter II. 

B. POISSON DEMAND CASE 

We shall now determine an exact expression for the 
TWUS^(s^) for the case in which demands are Poisson distributed. 
Let the mean rate of demand be A^ units per year and the 
lead time be a constant T^ . In addition to treating the 
demand variable as being discrete, the number of buys s^ 
also will be treated as a discrete variable. Thus if D^ is 
the lead time demand item j : 



n 



y c . s . < b 

i-1 1 1 - 



time weighted units short when there 
are s^ units for item i 



total sum of lead time demand 

( Jx X i T i‘ 

essentiality code for item i 
budget limit in a given period 



price of each item. 
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p(D i =s i } 



-X.T. s. 

i 1 1 ( X . T . ) 1 

l l 

(s i ) ! 



= p(s i ;X i T i ) 



( 1 ) 



Let 



p(s i ) = prob(D i >_s i ) 



d . =s . 

l l 



p(d i ;X i T i ) 



( 2 ) 



If there are units of stock for item i, Richards and 
McMasters [Ref. 8] show that the expected time-weighted 
units short in (0,T\) is given by 



E[TWUS i (s i ) ] 



T. 

-^•{p(s.+l) [X.T.-2s. + 

2 c i li l 



s. (s.+l) 

l i 

X .T. 
l l 






(3) 



For those cases where the expected lead time demand is 

large, the Poisson probabilities in (3) can be approximated 

by a normal distribution with mean X.T. and variance 
1 ii 

a? = X . T . . Let 

i ii 



<p (x) = -1— exp(-x 2 /2) 

/2tF 



be the standard normal probability density function and let 

00 

f(x) = / <f>(u)du be the complementary cumulative distribution 

x 
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function for the standard normal. Then expression (3) can 
be rewritten in terms of the normal probability function as 
follows : 



T. s.+l-A.T. 

E [ TWUS . ( s . ) ] = — 

ii 2 a . 



s . (s . +1) 

)|x iV 2s i + -- X -T — 

1 1 



+ T7 *' Ii a^ i >( X i T i- s i )) 



(4) 



This expression should be used in those cases in which A^T^ 
is large. We have developed the expression for the expected 
time-weighted units short in a period of length T^ when there 
are s. units of stock for item i. 

l 

In the next section we write the expression for the total 
essentiality-weighted time weighted units short over all items, 
and we provide a solution procedure for allocating the given 
budget optimally. 

C. SOLUTION PROCEDURE 

The mathematical program for the time-weighted-units short 
problem is: 

. N E.T. _ s . (s.+l) 

(Cl) min Z(s) = SIT I \.T. 

1=1 1 1 

+ pfs.^T.) (A.T i -s.)} 

N 

s.t. T c.s. < B 

i-i 1 1 - 
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To solve this problem we can use the Lagrangian multiplier 
technique. Let 



N N 

L(s. ,s~, . . . ,s ;0) = l Z . (s . ) + 9 (B - £ c.s.) (5) 

x n i=l 1 1 i=l 1 1 

Here Equation (5) is separable in the items, and minimization 
of the total objective function is accomplished by minimizing 
the individual functions Z^(s^) subject to budget constraints. 
Consider a single item i. Let 

AL i (s i ) = L^Sj-1) - L (s^) 

= z i (s i -i) - Z.( Si ) + 9c iS . - ec.(s.-l) 

= AZ^s.^) + Qc ± (6) 



where 

AZ . (s . ) = E. [TWUS. (s.-l) - TWUS . (s . ) ] (7) 

l l l ll li 

As shown earlier 

TWUS(s-l) = |{P(s) [AT-2(s-l)+ S( ^~ 1} ] +p ( s- 1 ; XT) (AT-s+1) } 

2 

= |{P(s) [AT-2s+^^-2J + ^p(s; XT) (AT-s+1) } 
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TWUS(s) = j{P(s+l) [AT-2s+Si|±iI + pCs;XT) (XT-s) } 



?{P(s) [AT-2s+ — * — ] -p (s ; AT) ( AT-2s+ - ( f^ 1} ) 



AT 



AT 



+ p ( s ; AT) (AT-s) } 



so that 



TWUS(s-l) - TWUS(s) 



£{P(S) [2 + s_ 2 - s - s2 - s 



2 

AT ] +p ( s ; AT) [s - frf + H r + AT_2s 



+ - AT + s] } 



= |[P(S) 12 “ ff] + f| p ( s ; AT) ] 



= P(s) [T - J) + j p (s ; AT) 



(8) 



Substitute Equations (7) and (8) into (6). Then 



s . _ s . 

AL . ( s . ) = E. [(T. - T ^)P(s. ;A.T.) + -^ p(s. ;A.T.) ] 

l i l i A^ ill A^ i i l 



+ 0c . 

l 



(9) 



The optimum solution s^ is the largest s^ such that 



AL . (s . ) > 0 

l i — 
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or equivalently, 





c . 

1 



s . 



+ ^7 P( s j_) ] 1 “ e 



The basic algorithm for solving this problem was explained 
in the previous chapter. A computer program for searching 
for 0 is provided in Appendix B. 

D. SAMPLE DATA RUNS 

Consider a weapon system which consists of 10 components. 
Suppose that the demand for each component is Poisson dis- 
tributed with parameter A ^ and the lead time is known constant 
T^. Let the budget available for procurement be $19224. 

Table III shows the optimal allocations provided by the TWUS 
model. 

The allocation given when the demand distribution is 
approximated by the normal distribution is also provided 
in Table III for comparison. (For comparability, the variance 
for the normal distribution is taken to be the same as the 
mean). Comparing the results, we observe that the normal 
case buys more of the high demands low cost items. There is 
a small difference in the allocation for items 8 and 9 which 
are more expensive than the others. For demand rates less 
than 10, the usefulness of the approximation is questionable. 
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TABLE III 



OPTIMAL ALLOCATION FOR TWUS MODEL 



Item 


X . 

l 


T. 

i 


c ± ($) 


E. 

i 


Poisson 


Normal 


1 


10 


1 


10 


1 


17 


20 


2 


100 


1 


20 


1 


113 


120 


3 


15 


1 


80 


1 


19 


21 


4 


20 


1 


2 


1 


32 


36 


5 


50 


1 


5 


1 


65 


71 


6 


80 


1 


30 


1 


90 


96 


7 


20 


1 


1 


1 


35 


37 


8 


15 


1 


200 


1 


17 


15 


9 


75 


1 


100 


1 


77 


74 


10 


10 


1 


75 


1 


14 


16 


The 


resulting 


# 

values of 


the 


ob j ective 


function for the 



optimal solutions are: 

Poisson Case Normal Case 



Z (s*) 

Shadow price 
Budget limit 
Budget left 



0.00094 
(6*) 0.00015 

$19224 
0 



0.0015 

0.00055 

$19224 

0 



The objective function value for the Poisson demand case is 
less than the normal demand case. The main reason for the 
difference is due to items 8 and 9 . 
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IV. PSEUDO AVAILABILITY MODEL 



A. DESCRIPTION OF MODEL 

In the previous chapter, the TWUS objective function was 
introduced as a means for allocating a limited budget. 
Operational availability is a widely stated measure of the 
operational readiness of military forces and weapon systems. 
Thus, it is appropriate to consider stockage models with an 
availability objective as a means of allocating limited 
resources . 

The most direct and meaningful measure of the influence 
of peace time operating stocks on readiness is weapon 
system (or end item) availability. We use the terms 
availability, end-item availability, and weapon system 
availability interchangeably to mean the probability 
that an end item, such as a tank or an aircraft, 
selected at random, is not waiting for a component 
to be repaired or shipped to it. [Ref. 6] 

Many authors have attempted to determine stockage levels for 

components by maximizing equipment operational availability, 

subject to a budget constraint. See, for example, Jee [Ref. 

7] . Usually, the availability for component i is defined by 

the ratio 




MTBF . 
1 



MTBF^ + MTTR i + MSRT^ (s^) 



where : 



MTBF. 



the Mean Time Between Failure of item i; 



l 
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MTTR^ = the Mean Time To Repair for item i; 
MSRT^(s^) = the Mean Supply Response Time for item i. 

If the weapon system is assumed to consist of the n components 
all arranged in series, then the system availability is the 
product of the individual item availabilities. (This assump- 
tion means that the system will fail if any of the components 
fails.) With this assumption, the allocation problem, stated 
in terms of system availability is: 



(PI) 



n 

max II A . ( s . ) 
s i=l 1 1 



n 

s.t. 7 c.s. < B i = 1,2,...,N 
i=l 



where : 



A. (s.) = the availability of item i having s. 

1 1 units in stock; 

B = the budget limit; 

= the price for each item i. 



In the 
cal of 
of the 
can be 



expression of A^(s^), the term MTBF^ is the recipro- 
the failure rate i, MTTR^ is assumed to be independent 
decision variables and the available funds and MSRT^(s^) 
expressed in terms of TWUS(s^) as 



MSRT i (s i ) 



A .T. 

l l 



TWUS.( S . ) . 
1 l 
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Thus, the main determination of availability from the point 
of view of the supply system is MSRT^(s^). Many techniques 
for solving this model have been developed. In the next 
section we represent an algorithm for solving the availa- 
bility model by using the marginal analysis method. 

B. SOLUTION PROCEDURE 

The model (Pi) is not additive in the individual com- 
ponent availabilities but is converted into an additive 
function by transforming the objective function. Taking the 
natural log of the objective function, the model can be ex- 
pressed in the following way. 

n 

(P2) max l in A . (s . ) 

si-1 11 

n 

s.t. 7 c.s. < B i = 1,2,. . . ,N 

. u , i l — 

1=1 

Now the model (P2) is separable for all i and maximization 
of (P2) yields the same solution as maximization of (Pi) . 

The marginal analysis method selects an item which gives at 
each step the greatest increase in log(A^(s^)) per dollar 
spent. 

STEP 1. Start with zero units for all items. 

STEP 2. Compute the increase in log availability per dollar 
spent as a result of purchasing one additional 
unit. 
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Jin (s^) ] 



A ^n A i ( s i ) 



c . 

1 




A. (s . ) - 

l l 



where i = 1 , 2 , . . . ,N 



Jin A . ( s . ) 

l l 



MTBF i 

£n t MTBF i + MTTR^ + MSRT^ (s^) ^ 



STEP 3. Select th&t item i corresponding to the maximum 
ratio . 



max 

all 



AJln A^(s^) AJln A 2 (s 2 ) 



AJln A (s ) 
n n 



c 

n 



STEP 4. Increase the number of units stocked for the item 
selected at step 3 by one additional unit if the 
unit price is less than the amount of budget 
remaining . 

STEP 5. Update the S vector, the MSRT(s) expression and 

decrement the available budget. If the remaining 
budget is greater than the cost of the cheapest 
item. Go to Step 3. Otherwise, Stop. 



In the following section, we will illustrate this procedure 
with a sample system. The computer program is provided in 
Appendix C. 
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C. A NUMERIC EXAMPLE FOR THE MSRT MODEL AND THE AVAILABILITY 

MODEL 

In the expression for availability, the MTBF and MTTR 
terms are not functions of the number of spare parts. There- 
fore it is commonly believed that maximization of system 
availability is equivalent to minimization of mean supply 
response time. However, this is not the case, as shown 
below. 

Suppose a weapon system consists of three components and 
the demands are Poisson distributed with parameters X^, 
and X^, respectively. The lead time is a known constant and 
the components have essentiality codes E^ . The unit price 
and MTTR are known and the budget is limited to 20 dollars. 
This information is summarized in Table IV. 

TABLE IV 

INPUT DATA FOR EXAMPLE 



ITEM 


A i 


C . 
1 


MTTR. 

1 


E i 


T. 

1 


1 


1 


5 


0.0274 


1 


1.0 


2 


0.1 


5 


0.0027 


3 


1.0 


3 


10 


1 


0.0054 


1 


1.0 



To solve MSRT minimization problems, we first determine 
the MSRT's for all possible cases. These values are provided 
in Table V. 
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TABLE V 



MSRT DATA FOR ALL FEASIBLE SOLUTIONS S. 

1 



MSRT ( s . ) 

l 


ITEM 1 


ITEM 2 


ITEM 3 


MSRT (0 ) 


0.9482 


0.9837 


0.6 


MSRT (1) 


0.3161 


0.1967 


0.5 


MSRT (2) 


0.0708 


0.02 


0.4099 


MSRT (3) 


0.0132 


0.00227 


0 .329 8 


MSRT (4) 


0.0021 


0.0002 


0.2596 


MSRT (5) 


0.0003 




0.1992 


MSRT (6) 






0.1485 


MSRT (7) 






0.1072 


MSRT (8) 






0.0746 


MSRT (9) 






0.0500 


MSRT (10) 






0.0322 


MSRT (11) 






0.0199 



Using the solution procedure described in the previous chap- 
ter we determine the optimal solution to be as shown in Table 
VI. 
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TABLE VI 



THE ALLOCATION OF SPARE PARTS FOR MSRT MODEL 



ITEM 



USED 



ALLOCATION 


AZ l (s l } 


AZ 2 (S 2> 


az 3 (s 3 ) 


BUDGET (?) 


(0,0,0) 


0.12462 


0.47216 


0.1 


0 


o 

O 


0.12462 


0.1040 


0.1 


5 


(1,1,0) 


0.04905 


0.1040 


0.1 


10 


(1,2,0) 


0.04905 


0.01265 


0.1 


15 


(1,2,1) 


0.04905 


0.01265 


0.09 


16 


(1,2,2) 


0.04905 


0.01265 


0.08012 


17 


(1,2,3) 


0.04905 


0.01265 


0.07022 


18 


(1,2,4) 


0.04905 


0.01265 


0.06039 


19 


(1,2,5) 


0.04905 


0.01265 


0.0507 


20 



The optimal solution for MSRT model is (1,2,5). Repeating 
the analysis for the availability objective function we 
obtain the results provided in Table VII from the marginal 
analysis procedure. 



AZ.(s.) = — [Z.(s.+1) - Z.(s.)] 

11 c^ 1 1 11 



— [in A . ( s . +1) - in A.(s.)] 

C . 11 11 

1 
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TABLE VII 



THE ALLOCATION OF SPARE PARTS FOR AVAILABILITY MODEL 
ITEM 123 

USED 

ALLOCATION AZ (s^ AZ 2 (s 2 ) AZ 3 (s 3 ) BUDGET 



(0,0,0) 


0.07712 


0.1294 


0.1529 




0 


o 

O 

** 

H* 


0.07712 


0.1294 


0.1611 




1 


(0,0,2) 


0.07712 


0.1294 


0.1689 




2 


o 

o 

u> 


0.07712 


0.1294 


0.1759 




3 


o 

o 


0.07712 


0.1294 


0.1808 




4 


(0,0,5) 


0.07712 


0.1294 


0.1821 




5 


VO 

o 

o 


0.07712 


0.1294 


0.1777 




6 


(0,0,7) 


0.07712 


0.1294 


0.1660 




7 


o 

** 

o 

00 


0.07712 


0.1294 


0.1469 




8 


(0,0,9) 


0.07712 


0.1294 


0.1217 




9 


(0,1,9) 


0.07712 


0.0447 


0.1217 




14 


(0,1,10) 


0.07712 


0.0447 


0.0938 




15 


(0,1,11) 

• 


0.07712 


0.0447 


0.06702 

• 




16 


• 

# 

(0,1,15) 






• 

• 




20 


The optimal 


solution for 


the availability model 


is 


(0,1,15) 


Comparing the results of 


the two models 


, we see 


that the 


availability 


model allocates more units 


to the 


high 


demand 


lower cost items than the 


MSRT model. 
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D. SAMPLE DATA RUNS 



Suppose that a weapon system consists of 10 components 
and the demand of each component is Poisson distributed with 
parameter A^, and lead time T^, mean time to repair MTTR^ 
are known constants. In order to maximize the availability 
of spare parts with budget constraint, we can use the modi- 
fied Availability model (P2) instead of (Pi) . By using the 
computer program in Appendix C, this problem can be solved. 
Table VIII provides the allocations of spare parts in the 
Availability model when the budget is 1170 dollars. 



TABLE VIII 

THE ALLOCATION OF SPARE PARTS FOR THE 
AVAILABILITY MODEL 



ITEM 


A . 

l 


T . 
l 


c i ($) 


E. 

1 


MTTR . 
1 


ALLOCATION 


A i (s i> 


1 


1.0 


1.0 


10.0 


1.0 


0.0137 


4.0 


0.986 


2 


0.1 


1.0 


20.0 


1.0 


0.0274 


2.0 


0.997 


3 


3.0 


1.0 


100.0 


1.0 


0.0137 


3.0 


0.821 


4 


25.0 


1.0 


2.0 


3.0 


0.0822 


37.0 


0.327 


5 


1.0 


1.0 


5.0 


1.0 


0.0274 


5.0 


0.973 


6 


0.5 


1.0 


5.0 


3.0 


0.0027 


4.0 


0.999 


7 


10.0 


1.0 


1.0 


1.0 


0.0054 


21.0 


0.949 


8 


5.0 


1.0 


100.0 


1.0 


0.0411 


3.0 


0.538 


9 


1.0 


1.0 


50.0 


1.0 


0.0082 


3.0 


0.987 


10 


2.0 


1.0 


100.0 


1.0 


0.1370 


2.0 


0.697 
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From the table, one can see that the availability of an 
item is greatly influenced by the MTTR term (see item 4) . 

The availability for that item never exceeds 0.333 even if 
the MSRT is zero. We also observe that the availability 
model tends to stock the high demand low cost items. 

The objective function for the optimal solution is given 

by : 

Total obejctive value 0.08999 

Shadow price 0.000261 

Budget limit $1170 

Budget left $0.0 

A comparison of the above results with the allocation given 
in Table VIII shows that the total availability is relatively 
low even though most of the items have high availabilities. 
Also as mentioned above, when the MTTR data for an item is 
large relative to the MTBF , a high availability cannot be 
achieved. 
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V. COMPARISON OF MODELS 



A. ANALYSIS FOR SAME DATA 

In this chapter, we continue to consider the allocation 
of spare parts to maximize the system performance in the 
different allocation models. In this thesis, we have looked 
at three models: the units short model, the time-weighted 

units short model, and the availability model. Since each 
model attempts to reduce stockouts as much as possible the 
allocations generated by the models are strongly correlated. 
This is especially true for the availability model and the 
MSRT model since availability is a function of MSRT. However 
we saw earlier that the allocations from the models are not 
necessarily the same. 

Assume that a weapon system consists of 10 items, the 

demands are Poisson distributed and MTTR. , c., T. , E. are 

1 1 i i 

known constants and a budget constraint of the weapon system 
is $1170. The optimal allocations for the three models are 
shown in Table IX. As can be seen, the TWUS model is more 
sensitive to the lead times than are the other two models 
(see items 5, 6, and 7). 

The units short model is more sensitive to the price of 
the item than are the other two models. For item 9 the units 
short model bought nothing, but the TWUS model and the availa 
bility model allocated 2 and 3 items respectively. All 
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TABLE IX 



THE ALLOCATIONS OF SPARE PARTS FOR THE 
THREE DIFFERENT MODELS 



Units 



Item 


X i 

(yr) 


cost 

($) 


Ess . 


T. 

i 

(yr) 


MTTR 

(yr) 


Short 

Model 


TWUS 

Model 


Avail . 
Model 


1 


1.0 


10 


1 


1 


0.0137 


3 


3 


4 


2 


0.1 


10 


1 


1 


0.0137 


2 


1 


2 


3 


15.0 


3 


1 


1 


0.0137 


24 


20 


22 


4 


15.0 


3 


3 


1 


0.0274 


26 


23 


24 


5 


3.0 


10 


1 


0.5 


0.0274 


4 


3 


4 


6 


3.0 


5 


3 


0.5 


0.0274 


6 


4 


6 


7 


10.0 


50 


1 


0.2 


0.0054 


0 


0 


3 


8 


10.0 


50 


1 


1 


0.0411 


8 


7 


2 


9 


2.0 


50 


1 


1 


0.0137 


0 


2 


3 


10 


2.0 


100 


4 


2 


0.1370 


5 


5 


5 



three models are highly affected by the essentiality code. 
This is illustrated by a comparison of items 9 and 10. Table 
X presents the corresponding values of the three objective 
functions . 
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TABLE X 



THE COMPARISON 


OF OBJECTIVE 


VALUES FOR 


THREE MODELS 


OBJ 








FN 


UNITS 






MODEL 


SHORT 


TWUS 


AVAILABILITY 


UNITS SHORT 


0.1523 


0.0458 


0.0483 


TWUS 


0.1527 


0.0353 


0.0687 


AVAILABILITY 


0.1969 


0.0775 


0.0709 



The above table was established by computing each objective 
function for the allocations determined by the three differ- 
ent procedures. Comparing the results of the three models, 
the TWUS model seems to do the best job considering all three 
objective functions. However, no general conclusions can be 
drawn about the preference of the TWUS model for other 
situations . 

One needs to determine which objective function most 
closely matches a servicers' feeling about how operational 
readiness is affected by stockouts and delays in satisfying 
stdckouts . 

B. DISCUSSION OF SIMILARITIES 

In the budget allocation problem there are many factors 
which affect the allocation such as demand, lead time, cost, 
time to repair and essentiality. The three models share simi- 
lar properties. First of all, as can be seen in the above 
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example, all models tend to stock the cheap, high demand 
items in favor of expensive low demands items. This is 
because of the models attempts to get the biggest benefit 
per dollar spent. Potential benefit per additional unit 
increases with an items demand rate. Second, items having 
high essentiality code are given preference, as is the in- 
tent of essentiality assignment schemes. Essentiality 
weighting is one way to counter the preference given the high 
demand low cost item observed earlier. It is frequently the 
case that the most critical items are low demand expensive 
items. Without the essentiality weighting such items would 
be neglected by the type of models examined in this thesis. 
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VI. USE OF THE MODELS FOR BUDGET DETERMINATION 



A. EFFECTIVENESS VS. BUDGET 

The models that we have discussed have attempted to opti- 
mize performance subject to a budget constraint. We have 
assumed that the budget was given. There are many ways in 
which budgets are determined. However, budgeting people and 
inventory managers alike often express the desire to have a 
methodology that they can use to determine the amount of 
money that should be provided. 

In most cases the amount is determined historically by 
giving an amount equal to what has been provided in the past 
for similar systems or perhaps by giving a little more or 
less based on judgement or financial constraints. There is, 
however, a strong interest brought about by Congressional 
pressures to relate resources to readiness. Congress wants 
to know "how much money is needed to support our weapon sys- 
tems at a specified level of performance." In this chapter we 
show how the models developed earlier in this thesis can be 
used in just this manner. 

Specifically, we show how the models that we have developed 
can be modified easily to determine the minimum budget re- 
quired to provide a specified level of logistics performance. 

The models developed earlier can each be run for a range 
of budget levels producing for each given budget an allocation 
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and a predicted overall level of performance. Figure 3 illus- 
trates this for the case in which the performance measure is 
pseudo-availability. As expected, the curve shows that 
availability is a non-decreasing function of budget with 
decreasing marginal returns. This can be done also for the 
time-weighted units short model or any of the other models 
discussed in this thesis. In all cases we would obtain a 
similar display. Performance is a monotonic function of 
budget with decreasing marginal returns. 

Figure 4 displays a similar result for the case in which 
the performance measure is MSRT. Each point on the curve 
represents an optimal level of performance for a given budget. 
For this example displayed, Figure 3, there is a little bene- 
fit to be gained by increasing the budget above $2500. How- 
ever there is a dramatic increase in effectiveness obtained 
by increasing the budget from $1000 to $2000. This is pre- 
cisely the sort of information needed to make intelligent 
budgeting decisions. Of course some decision maker must 
decide if the increase in effectiveness is worth the addi- 
tional expenditure. 

If a specific level of effectiveness is specified, one 
can graphically determine the amount of budget required by 
simply moving horizontally across the graph from the speci- 
fied level of effectiveness until the curve is intersected 
and then down to the budget axis. 
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AVAILABILITY VS 0UOGFT 




LEGENO 
o * X 



Figure 3. Total Availability Vs. Budget Curve 
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MSRT VS BUDGET 




LEGEND 

OaWC 



Figure 4 . MSRT Vs . Budget Curve 
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The next section determines analytically the minimum 
amount budget required by solving a companion problem to 
the problems discussed earlier in this thesis. 

B. COMPANION PROBLEMS 

In the previous chapters we have concentrated on the 
optimization of system effectiveness with a budget constraint. 
For many weapon systems such as air detection radars, missiles 
and nuclear delivery systems, the system performance is so 
important that the necessary budget will be provided to 
attain whatever performance is deemed necessary. 

For such systems it is reasonable to restate the optimi- 
zation problem to determine the minimum budget required to 
satisfy a specified level of performance. Consider, for 
example, the availability optimization problem and the com- 
panion problem: 

n 

(Dl) min J c.s. 

i=l 1 1 

n 

s.t. II A. ^ L i = 1 ,2 , . • . ,n 

i=l 1 

where L is the minimum performance level for a weapon system. 
For the MSRT model case the corresponding problem is: 

(D2) min £ c i s i 

n 

s.t. T MS RT . ( s . ) < R i = 1 , 2 , . . . ,n 

i-i ii- 
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where R is the maximum allowable cumulative supply response 
time for the weapon system. 

Problems (Dl) and (D2) can be solved using the same 
methods explained in Chapters III and IV. In the above 
models the budget is determined so that the system require- 
ment for availability or main supply response time can be 
achieved. 

For problem (D2) the total cost is minimized when 
n 

7 MSRT. (s . ) = R . 

. L , 11 

i=l 

Sometimes a minimum allowable supply response time is 
required for each item. In such a case multiple constraints 
could be specified. This is illustrated below: 

n 

(D3) min Y c.s. 

i=l 1 1 

MSRT^s^ < R 1 

msrt 2 (s 2 ) <_ r 2 
# 

MSRT (s ) < R 

n n — n 

where R^ is the maximum allowable supply response time, 
i = 1,2,. . . ,N. 
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To solve (D3) , find the smallest s^ such that 



MSRT 1 (s 1 ) 



R 



1 



msrt 2 (s 2 ) 



R 



2 



MSRT ( s ) 
n n 



R 



n 



This problem is solved easily using the same procedures 
which we discussed. 

So far we have discussed many different ways to apply 
the theoretical models to practical use of models for budget 
determination decisions. There is no unique method which 
gives us an optimal result. So the user of these models 
should choose one of possible methods so as to maximize the 
system performance or minimize the total cost. 
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VII. CONCLUSIONS 



It is concluded that the various measures of effectiveness 
can be used in the budget constrained multi-item inventory 
system with stochastic demands. We have examined some of 
the more reasonable measures like minimization of units short, 
minimization of time-weighted units short and maximization 
of system availability. We have also looked at models which 
incorporate essentiality weights into each of the models. 

In order to solve budget allocation stockage problems a 
feasible, efficient method of effecting line item inventory 
control is available using an adaptation of Everett's 
Generalized Lagrangian Multiplier method. Further, the use of 
a G.L.M. procedure provides valuable information for system 
managers as to relative effectiveness of additional procure- 
ment funds, versus additional transition processing capability. 
The final value of Lagrangian multipliers can be interpreted 
as the amount of improvement of the objective function per 
unit dollar spent. 

The models discussed in this thesis are all more likely 
to stock cheap, high demand items than expensive, low demand 
items. Such is the nature of budget constrained optimization 
problems. If a system manager wishes to maintain enough 
stock for an item having low demand, high cost, his only 
alternative in our models is to assign a high essentiality 
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code for the item. The essentiality code has the effect of 
reducing the ratio C/E as opposed to C. In the solution 
procedure for each model, the assigned essentiality code 
directly affects the allocation for the item. 

We have shown how the models can be used as a tool to 
determine the amount of budget. A simple graphical procedure 
allows a decision maker to determine the minimum budget re- 
quired to search a specified level of performance. The opti- 
mization model is run several times to generate a plot of 
performance vs. budget. Each point of the curve represents 
the effectiveness for the optimal allocation of a given budget. 
A manager can, first, determine an appropriate system per- 
formance level and read from the curve the budget required 
to achieve the effectiveness. 

Further analysis to improve these models may be possible. 
For instance, it would be useful to have an automatic search 
algorithm for the Lagrangian multipliers for a multiple 
constrained problem. It may also be possible to relax the 
assumptions for a constant lead time or a constant mean time 
to repair. These single period inventory may expand to time- 
dependent multi-item, multi-echelon, multi-indenture inventory 
systems . 
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COMPUTER PROGRAM FOR INTERACTIVE SEARCH 
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APPENDIX B 



COMPUTER PROGRAM FOR AUTOMATING SEARCH 
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